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Abstract 

HF, H 2 0, CN - and their hydrogen-bonded complexes have been studied using 
state-of-the-art ab initio quantum mechanical methods. A large Gaussian one- 
particle basis set consisting of triple zeta plus double polarization plus diffuse s and 
p functions (TZ2P Hr diffuse) was used. The theoretical methods employed include 
self-consistent-field, second-order Moller-Plesset perturbation theory, singles and 
doubles configuration interaction theory and the singles and doubles coupled cluster 
approach. The FH~CN", FH* ~NC- and H 2 0"CN-, H 2 0 ”NC- pairs of complexes 
are found to be essentially isoenergetic. The first pair of complexes are predicted to 
be bound by ~ 24 kcal/mole and the latter pair bound by - 15 kcal/mole. The ab 
initio binding energies are in good agreement with the experimental values. The two 
pairs of complexes exhibit small structural differences with the N"’H hydrogen bond 
being shorter than the analogous C~H hydrogen bond. The infrared (IR) spectra 
of the two pairs of complexes are also very similar, though a severe perturbation of 
the potential energy surface by proton exchange means that the accurate prediction 
of the band center of the most intense IR mode requires a high level of electronic 
structure theory as well as a complete treatment of anharmonic effects. The bonding 
of anionic hydrogen-bonded complexes is discussed and contrasted with that of 
neutral hydrogen-bonded complexes. 
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Introduction 


Over the past sixty years, hydrogen-bonded complexes have attracted consid- 
erable attention from chemists. Much of the interest has been directed at the under- 
standing of the nature of the relatively weak bonding present in neutral hydrogen- 
bonded complexes. To this end, several different hydrogen bonding decomposition 
schemes have been developed. The basis for the classical description of hydrogen 
bonding was presented in a review 1 by Coulson in 1957. The classical hydrogen 
bond energy is decomposed into four distinct components - (1) the electrostatic 
energy; (2) the delocalization energy (commonly referred to as the charge transfer 
energy); (3) the repulsive energy and (4) the dispersion energy. Since Coulson lim- 
ited his review to hydrogen-bonded complexes involving a polar molecule containing 
an electronegative atom (such as N, 0 or F) and a molecule containing a polar A-H 
bond (where A = N, 0 or F), the electrostatic interaction is viewed as the dominant 
attractive force. For some Van der Waals complexes Morokuma and coworkers have 
demonstrated 2-4 that the electrostatic energy may be very small or even represent 
a repulsive force. However, for most hydrogen-bonded complexes the electrostatic 
interaction will be attractive. The dispersion energy also represents an attractive 
force and thus, in Coulson's review, the “repulsive force” is the only interaction 
which separates monomers A and B. The explanation of the physical nature of 
this “repulsive force” is based, not surprisingly, on electron-electron repulsion, i.e. 
the mutual repulsion of the electron cloud of monomers A and B, and quantum 
mechanical effects are not discussed. 

Subsequently Morokuma and coworkers 2 4 extended and adapted this decom- 
position scheme into a rigorous quantum mechanical approach as viewed through 
the self-consistent-field (SCF) ab initio method. There are six components in this 
decomposition scheme - 1) electrostatic; 2) polarization; 3) exchange repulsion; 4) 
charge transfer; 5) “MIX”; and 6) “CORR.” The CORR term is the contribution 
of electron correlation which Morokuma and coworkers did not investigate in detail 
though they stated that the most significant portion of the intermolecular correla- 
tion energy is known as the dispersion energy which is an instantaneous effect due 
to the simultaneous correlation of electrons in monomer A and monomer B. The 
MIX term is the higher order couplings of the first four components. The polariza- 
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tion interaction is the distortion of the electron density of A (B) by the presence 
of monomer B (A) and higher order effects. In applying this decomposition scheme 
to normal hydrogen-bonded complexes, Umeyama and Morokuma 4 concluded that 
the binding in these chemical systems is mostly electrostatic in nature with a small 
but significant contribution from the charge transfer energy. 

More recently qualitative approaches based upon electrostatic and polariza- 
tion interactions have been developed for the theoretical prediction of molecu- 
lar structures 5 ’ 6 and vibrational frequency shifts 7 of hydrogen-bonded complexes. 
When applied to neutral hydrogen-bonded complexes, both of these methods yield 
qualitatively correct results, although their accuracy is generally not quantitative 
and, in some cases, not even semi-quantitative 5 ' 7 . Furthermore, based upon the 
results of these approaches which have been reported thus far, it seems likely that 
these classical, perturbative approaches will break down as the binding energy of 
the complex increases. Since anionic hydrogen-bonded complexes are typically much 
more strongly bound, these simple approaches are not likely to be as successful. 

A more rigorous approach to the study of weakly bound systems which has 
been applied with much success has been the use of ab initio quantum mechanical 
methods 5 . Numerous studies have demonstrated that the SCF method (coupled 
with a large one-particle basis set) is capable of describing a weak hydrogen bond 
reasonably well, except for the dispersion energy. However, the recent formulation 
and development of better and more efficient electronic structure methods has en- 
abled the direct quantum mechanical investigation of weakly bound molecular com- 
plexes at correlated levels of theory. For example, Handy and coworkers 9 ' 12 have de- 
termined the equilibrium structures of several weakly bound hydrogen-bonded com- 
plexes [HCN-HF, HCN-HC1, (C 2 H 2 ) 2 , (C 2 H 2 ) 3 , FH-CO and FH ’ NNO] using 

large one-particle basis sets in conjunction with second-order Mpller-Plesset pertur- 
bation theory (MP2) and have found good agreement with experiment. However, as 
Rice, Lee and Handy (RLH) have demonstrated 12 with their study of H 2 CO "HC1, 
MP2 is not always adequate, especially when electron correlation effects are very 
important in the binding of the complex. RLH found that the theoretically more 
complete coupled-pair functional (CPF) approach gives much better H 2 CO and HC1 
monomer properties and, consequently, the H a CO~HCl structure is in excellent 
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agreement with the limited experimental data. In particular the dipole moment 
of HC1 is much better described with the CPF approach, supporting the thesis 
that electrostatic interactions are important in hydrogen bonding. Nonetheless, the 
substantial differences in the equilibrium geometry of the H 2 CO'”HCl complex ob- 
tained at the SCF, MP2 and CPF levels of theory demonstrate the importance of 
the dispersion energy. 

Theoretical and experimental studies of anionic hydrogen-bonded complexes 
axe more recent, especially in the gas phase. With the aid of three theoretical 
studies 13-15 , Kawaguchi and Hirota 16 have recently detected and analyzed the first 
high resolution infrared (IR) band of an anionic hydrogen-bonded complex (FHF ). 
There have been several theoretical studies of anionic hydrogen-bonded complexes, 
though very few of these have determined equilibrium structures and molecular 
properties beyond the SCF level of theory. Furthermore, none of the theoretical 
investigations have studied the decomposition of the hydrogen bond energy of an 
asymmetric anionic hydrogen-bonded complex. Umeyama et al. have performed 1 
a decomposition of the hydrogen bond energy of FHF - and find, not surprisingly, 
that charge transfer is much more important than for neutral hydrogen-bonded 
complexes. However, the decomposition analysis of FHF - is almost certainly not 
representative of asymmetric systems since FHF - adopts a D x h equilibrium struc- 
ture. 

In some respects anionic hydrogen-bonded complexes provide more of a chal- 
lenge than neutral and cationic hydrogen-bonded complexes for both experimental- 
ists and theoreticians. For example, the high resolution IR spectroscopist must deal 
with the very small population of anions that can be generated. Moreover, once 
a sufficient population has been attained, the analysis of the spectrum is further 
complicated by the presence of many other ionic species. The difficulty in the ab 
initio study of anionic species is well documented (see for example references 18-23). 
This difficulty generally arises due to the greater importance of electron correlation 
in anionic species. 

However, in other respects the study of anionic complexes is much easier than 
the study of similar cationic complexes. From an experimental viewpoint, the large 
binding energies of anionic hydrogen-bonded complexes should make their gener- 
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ation an easier task. In order to assess the implications of theoretical studies of 
anionic complexes, consider H 5 0+ and H 3 0 2 . These two systems are isoelectronic 
and it is likely that electron correlation effects will be more important for H 3 0 2 
than for H 5 0^". However, H 3 0^ has two fewer nuclei and therefore six fewer nu- 
clear degrees of freedom. The significant point is that most of the nuclear degrees 
of freedom which have been eliminated are large amplitude motions, an adequate 
treatment of which requires knowledge of a large portion of the potential energy 
surface (PES) as well as a sophisticated treatment of the nuclear motion problem. 
Existing methods for the accurate determination of the vibrational energy levels of 
polyatomic species which go beyond the harmonic oscillator approximation and are 
capable of adequately treating large amplitude motions are highly dependent upon 
the number of large amplitude nuclear degrees of freedom. Thus, while the de- 
scription of the electronic structure of anionic hydrogen-bonded complexes is more 
difficult, the accurate solution of the nuclear motion problem should be more feasi- 
ble. Therefore, the results of the current study provide data which will ultimately 
enable the detailed theoretical investigation of molecular systems with several large 
amplitude nuclear degrees of freedom. 

For anionic hydrogen-bonded complexes, the most difficult region of the PES 
to describe theoretically is the proton transfer coordinate which corresponds to the 
process AH + B~ -* A - + HB. The difficulty arises due to the possible exis- 
tence of two minima corresponding to A - 'HB and AH 'B . Several studies have 
investigated this region of the PES for symmetric and asymmetric anionic hydrogen- 
bonded complexes (see for example references 24-27). In addition, one of these*' 
also examined the adequacy of various vibrational analysis techniques. These stud- 
ies have demonstrated that electron correlation effects are vitally important 27 in 
obtaining a reliable description of the PES along the proton transfer coordinate, and 
that when A and B are both very electronegative atoms there is generally no barrier 
(and hence no second minimum) to proton transfer for asymmetric systems 

To date, no high level theoretical investigations of anionic hydrogen-bonded 
complexes involving HF and CN" or H 2 0 and CN~ have been reported. Experi- 
mentally, the IR spectrum of the M + FHCN" ion pair (M + being an alkali metal 
cation) has been studied via matrix isolation techniques by Ault . Fundamental 
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vibrations were observed in the 1100 cm 1 , 1800 cm 1 and 2500 cm 1 regions 
and assigned to a bending mode, the proton transfer mode (mostly H-F stretch) 
and the C-N stretching mode, respectively. Ault also observed that the 1100 cm -1 
band split into two components which were attributed to the presence of the metal 
cation M + . All three modes varied somewhat depending upon the composition 
of the matrix and the reactants used to form the M + FHCN ion pair. We note 
that Ault does not seem to have considered the existence of the FHNC isomer. 
Larson, McMahon and Szulejko 24 ’ 29 have determined the binding energies (i.e. hy- 
drogen bond strength) of both the FH - "CN - and the H 2 0 CN - complexes. The 
former is more strongly bound (21.1 kcal/mole) while the latter’s binding energy 
(12.7 kcal/mole) is still much larger than that of a typical neutral hydrogen-bonded 
complex. 

The purpose of this study is to obtain a better understanding of the electron 
correlation requirements in the ab initio study of asymmetric anionic hydrogen- 
bonded complexes and a more complete understanding of the nature of the bonding 
present in such systems. Thus, the conclusions of the present study will be useful 
in deciding upon the level of ab initio theory necessary to determine accurately the 
PES of an anionic hydrogen-bonded complex. The theoretical approach is described 
in the next section. The following sections contain a presentation and evaluation of 
our results. Concluding remarks are presented in the final section. 

Theoretical Approach 

It is well known 18-21 that large basis sets are necessary in order to obtain 
highly accurate results for anionic systems. Therefore, a single, large one-particle 
basis set has been used in this study. This basis consists of Dunning s (5s3p) 
contraction of Huzinaga’s 31 [10s6p] Gaussian primitive set for the heavy atoms (C, 

N, O and F). For hydrogen, the standard (3s) contraction 30 of the [5sj primitive 
set 31 was used. The hydrogen s function exponents were scaled by a factor of 1.49, 
as suggested by Dunning. In order to describe better the anionic nature of these 
systems diffuse s and p functions were added to the heavy atom basis (a 4(P (C) = 

O. 04812, 0.03389; a, tP (N) = 0.06742, 0.04959; a., p (0) = 0.08993, 0.05840; a, )P (F) = 
0.1164, 0.07161) while a diffuse s function was included in the hydrogen atomic basis 
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(a,(H) = 0.06696). These orbital exponents were determined in an even tempered 
manner using a method suggested previously 18 . Finally, two sets of polarization 
functions were added to all the atomic basis sets. The orbital exponents of the d 
polarization functions are a d = 1.5, 0.35 for the heavy atoms and a p = 1.4, 0.25 for 
hydrogen. These are the values suggested by van Duijneveldt 32 and have been used 
previously 18 in the study of anionic systems. This basis set is designated TZ2P -f 
diffuse. In all cases, the full complement of six Cartesian d functions was included 
in the basis giving 110 basis functions for the two larger complexes and 100 basis 
functions for the two smaller complexes. Linear dependency tests of the one-particle 
basis set were performed routinely and no problems were encountered. 

The first ab initio method utilized is the simplest, namely the restricted 
Hartree-Fock (RHF) SCF technique. As discussed previously, electron correlation 
effects must be included in order to account for the dispersion energy in hydrogen 
bonded complexes. Moreover, it is reasonable to expect that the dispersion energy 
(or electron correlation effects) will be more important for anionic hydrogen-bonded 
complexes because of the diffuse, polarizable nature of the electron cloud of anions. 
Therefore, three different electron correlation methods have been used in order 
to investigate the importance of electron correlation effects. The first approach 
is second-order Mpller-Plesset perturbation theory 33 (MP2). The second, another 
commonly used method, is singles and doubles configuration interaction (C1SD) 
which is based upon the variational principle. The third, and theoretically most 
complete method, is the singles and doubles coupled cluster approach (CCSD). The 
MP2 and CCSD methods have the advantage of being exactly size extensive and size 
consistent 34 . The CISD technique includes configuration mixing which MP2 does 
not take into account, but CISD is an n 6 procedure (MP2 is an n 5 procedure) where 
n is the number of active molecular orbitals. The CCSD method does allow config- 
uration mixing but is somewhat more expensive than CISD 22 . Also, the energetics 
of complicated chemical reactions are more easily computed using size consistent 
methods since “super-molecule” energies are not necessary. Therefore, based upon 
the above discussion and previous results 35 the CCSD method is expected to yield 
the most reliable results. 

Equilibrium structures of the complexes have been obtained with the SCF, 


7 


MP2 and CISD methods. Due to the computational cost and the available compu- 
tational facilities, at the CCSD level of theory it was only possible to optimize the 
FH-CN" and FH-NC - complexes. However, single point CCSD energies at the 
MP2 and CISD equilibrium structures have been performed for the H 2 0-CN" and 
H 2 0-NC~ pair of dimers. Also, in order to reduce the CISD expansions the heavy 
atom ls-like core molecular orbitals were required to be doubly occupied in all 
configurations and the corresponding virtual counterpart was deleted from the pro- 
cedure. The same procedure was also used in the CCSD optimization of the two 
FH;CN" complexes. 

As noted in the introduction, the structures of many neutral hydrogen-bonded 
complexes are strongly dependent upon the respective monomer properties. Thus, 
in an attempt to judge better the reliability of the theoretical predictions of the 
complexes, all possible monomers have been studied using the basis set and ab 
initio methods described above. 

In most cases, analytic energy gradient methods 36 have been employed to lo- 
cate precisely the equilibrium structures. Analytic energy second derivative meth- 
ods have been used to determine the SCF 37 and MP2 38 ’ 39 Hessian matrices while the 
CISD and CCSD Hessians were obtained numerically by taking central differences 
of analytic gradients. Infrared intensities have been determined via the double har- 
monic approximation. The dipole derivatives were determined analytically at the 
SCF level of theory and central differences of dipole moments were utilized at the 
CISD and CCSD levels of theory. In all cases, dipole moments were determined with 
respect to the center of mass and evaluated as energy derivatives 40 . In the numer- 
ical central difference procedures, energy invariance relationships for the Hessian 41 
and dipole derivative 42 matrices were used in order to reduce the number of gradi- 
ent evaluations. This is the most efficient numerical procedure for the evaluation of 
dipole derivatives provided that the numerical Hessian is also required. 

All SCF and MP2 investigations were performed with the Cambridge Analytic 
Derivatives Package 43 (CADPAC), while the CISD 44 ' 45 and CCSD 46 ’ 47 studies were 
performed with the Berkeley suite of programs modified to run on a Cray X-MP. 
The CCSD studies of the dimers were performed with a recently developed 22 vector- 
ized CCSD method. SCF, MP2, CISD and monomer CCSD calculations were per- 


8 


formed at the University of Cambridge. The CCSD optimizations of FH’CN and 
FH "NC - were performed on the Cray X-MP/48 at NASA Ames Research Center. 


Monomer Properties 

The equilibrium structures, total energies and dipole moments of the various 
monomer fragments are presented in Table 1. The CCSD method provides much 
better agreement with experimental structures and dipole moments. In fact, for 
the neutral molecules the magnitude of the errors in the CCSD prediction of the 
equilibrium structures are less than half those present for the CISD and MP2 struc- 
tures. The only exception arises for the bond angle in H 2 O where the MP2 result is 
fortuitously in better agreement with the experimental value. However, the CCSD 
value is only 0.3! too large. Nonetheless, the main conclusions to be drawn from 
the results of Table 1 are that the CCSD method, as expected, performs better 
than either the CISD or MP2 approaches and, more importantly, that for many 
chemical systems quantitatively accurate structures (i.e., Ar e < 0.001 Aand A0 e 
< 0.5° ) may be obtained with the CCSD electron correlation procedure coupled 
with a large one-particle basis set. 

Another important aspect concerns the equilibrium structure of OH . Note 
that the absolute magnitude of the error of the CCSD bond length is significantly 
larger than for the A-H (A = C,0,N) bonds of the neutral molecules. The electron 
correlation energy for anions is generally larger than for isoelectronic neutrals (as 
evidenced here by comparing the correlation energies of OH and H 2 O). Thus, 
anions usually require a more rigorous treatment of electron correlation in order to 
obtain accuracy comparable to that obtained with neutral molecules. With these 
considerations, it is not too surprising that the CCSD bond length of OH is not 
as accurate (compared to experiment) as the CCSD O-H bond length in H 2 O. 

Based upon previous experience, the only other geometrical parameter which 
is potentially difficult for ah initio methods is the C-N triple bond present in CN , 
HCN and HNC. The results given in Table 1 confirm the inherent difficulty in 
adequately treating the C-N triple bond, though again the CCSD equilibrium values 
for HCN and HNC are superior to either the CISD or MP2 quantities. For these 
three molecules, the C-N bond distance decreases in the order CN > HNC > HCN, 
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which yields insight into the nature of the carbon and nitrogen lone pairs. Since the 
C-N bond distance in CN - is longer than that in C-N radical, the HOMO exhibits 
slightly anti-bonding character. Thus, when a proton is attached to form either 
HCN or HNC the C-N bond distance shrinks due to the polarization of electron 
density away from the C-N linkage. Therefore, since the C-N distance in HCN is 
shorter than that found in HNC, we may conclude that C contributes more to the 
antibonding characteristics than does N. As an aside to the above discussion it is 
interesting to note that the CISD and CCSD correlation energies for these three 
molecules increase in magnitude in the order HNC < CN - < HCN. At the MP2 
level of theory the correlation energy of CN is slightly larger than that for HCN. 

The dipole moments of the monomers are predicted to almost equal accuracy 
with the CISD or CCSD methods, though the CCSD dipole moment is usually 
in better agreement with experiment. A noteworthy point which has particular 
relevance to this study is the fact that the rather sizable dipole moment of CN - 
(0.64 D) has carbon at the negative end. Consequently, there are competing effects 
as to which lone pair of electrons (the C lone pair or the N lone pair) will act as 
the better Lewis base, or in other words, which end of CN will form the stronger 
hydrogen bond? The nature of these phenomena may be understood by considering 
electron density maps of the C and N lone pair molecular orbitals which have been 
given by Taylor et al The C lone pair orbital is broad and diffuse whereas the 
N lone pair orbital is tighter. Thus, the C lone pair electrons will produce a larger 
attraction on the proton of the hydrogen-containing monomer in the AH'CN 
complex whereas for the AH'"NC — complex the N nucleus will have a stronger 
interaction with the electron cloud of the hydrogen-containing monomer since it 
will be able to approach more closely (i.e. form a shorter hydrogen bond). Thus, it 
is not evident, a priori , which of the two complexes will be more stable. Therefore, 
if CN - is one of the monomers of a hydrogen-bonded dimer, it will be necessary to 
investigate both AH 'CN - and AH 'NC - . As we shall demonstrate, for AH being 
either HF or H 2 O both sets of isomers are nearly isoenergetic. 

The harmonic vibrational frequencies and infrared (IR) intensities of the 
monomers are reported in Table 2. As has been noted by several authors recently, 
the CISD method seems incapable of properly describing the curvature of the PES 
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around an equilibrium point whereas the size extensive MP2 and CCSD methods 
both yield quite good harmonic vibrational frequencies. The CISD harmonic fre- 
quencies are consistently too high even with the rather large one-particle basis set 
used in this study. This particular inadequacy with CISD is believed to be related 
to the lack of size extensivity 49 , though no direct proof has as yet been given. 

Somewhat surprisingly the MP2 and CCSD harmonic frequencies are about 
equally accurate for this set of molecules (with the TZ2P+diffuse basis set), with 
the CCSD harmonic frequencies for OH", HCN and HNC being somewhat better 
than the MP2 values and the CCSD harmonic frequencies for HF and H 2 0 being 
marginally worse than the MP2 quantities. In any case, the MP2 and CCSD har- 
monic frequencies for the monomers are in very good agreement with experimental 
values with the possible exception of the bending mode in HCN and HNC. This 
particular normal mode is very sensitive to specific basis set deficiencies and the 
interested reader is referred to references 35 and 50 for more details of this effect. 
Based upon the CCSD and experimental results of the C-N stretching normal mode 
of HCN and HNC, the experimental harmonic frequency of CN~ can be estimated 
to lie near 2076 cm -1 . Using the u> e x e = H-3 cm -1 determined by Taylor and 
coworkers 48 , the experimentally unknown fundamental frequency is predicted to lie 
at 2053 cm -1 . This value is in excellent agreement with the high level calculations 
of Botschwina 51 (2052 ± 6 cm -1 ). 

The C-N stretch harmonic frequency decreases in the order HCN > CN" 
> HNC. Based upon the previously discussed C-N bond distances, the CN har- 
monic frequency would probably have been expected to be the lowest. This result 
demonstrates that caution must be exercised in relating geometric and vibrational 
properties. 

The IR intensities reported in Table 2 are consistent 52 with the expectation 
that electron correlation tends to reduce the magnitudes. The CCSD IR intensi- 
ties demonstrate that while CISD IR intensities are a vast improvement over SCF 
quantities, the CISD procedure still underestimates the correlation contribution to 
IR intensities. This observation is entirely consistent with a recent study 49 on the 
effects of triple and quadruple excitations in the Cl electron correlation procedure 
where it was shown that, like the electronic energy, many molecular properties 
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tend to converge from one direction as the excitation level is increased (i.e., do not 
exhibit oscillatory convergence). Thus, based upon the IR intensity and dipole mo- 
ment data, we may conclude that the CCSD approach better describes the electrical 
properties of the molecular systems included in this study. 


Energetics 

In order to determine which dimer will represent the global energy minimum 
(e.g. FH "CN - or F“"HCN), it is necessary to consider the enthalpy of the two 

reactions 

HF + CN" + AH 298 -> F" + HCN (1), 

and 

H 2 0 + CN - + AH 298 -* 0H~ + HCN (2). 

If AH 298 is positive, then the dimer will correspond to the reactants. Table 3 
contains ab initio and experimentally derived values for AH 29 8- The data in Table 
3 clearly indicate that the reactants of equations (1) and (2) should form the more 
stable dimer. This situation arises due to the large electron affinity of CN and the 
large F-H and 0-H bond energies. Therefore, as discussed above, we must consider 
two sets of isomers corresponding to hydrogen bond formation through the C or N 
end of CN - . 

Table 4 contains the total energy, binding energy and dipole moment of each 
complex, determined at the equilibrium structure for the given level of theory. The 
most important point to notice is that both sets of isomers are nearly isoenergetic. 
Thus, it is not possible to say definitively which stationary point represents the 
lowest energy structure. However, the ab initio data are reliable enough to con- 
clude that the actual difference between the isomer’s binding energies will not be 
greater than 5 kcal/mole. Thus, depending upon how the complexes are formed, 
it is possible that both isomers will be present under a given set of experimental 
conditions. 

The binding energies in Table 4 have incorporated a correction for the basis 
set superposition error (BSSE). The BSSE was determined using the counterpoise 
method at the SCF and CCSD levels of theory . It is well established 8 that the BSSE 
is generally larger at a correlated level of theory and that in order to reduce the 
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BSSE at a correlated level of theory a very large one-particle basis set is required 8 ’ 10 
. Therefore, we include the CCSD BSSE for the energetics determined with electron 
correlation methods. 

The best theoretical estimates of the [H20”CN - ; H20'"NC“] binding en- 
ergy (14.5 and 14.7 kcal/mole, respectively) are in good agreement with the 
experimental 29 value, 12.7 ± 0.8 kcal/mole. The agreement between theory and 
experiment for the binding energy of the [FH’"CN~; FH"NC~] set of isomers is 
also good, again being somewhat too large. However, the difference between the- 
ory and experiment is somewhat larger for the FH;CN~ pair of complexes. The 
experimental value may be somewhat too low for this complex and support for this 
assertion is found by comparing the theoretical and experimental binding energies 
of the F - “H 2 0 complex. A similar 27 level of theory to that used in this study gave 
a binding energy of 23.2 kcal/mole for F _ "H 2 O with the experimental quantity 
being 23.3 kcal/mole. Thus, the results of this study suggest that the FH;CN _ pair 
of complexes may be slightly more strongly bound than F~'"H 2 0, but experimental 
values suggest the opposite situation. As we shall show, other molecular properties 
determined via ab initio methods (such as the IR intensity of the proton trans- 
fer mode) are consistent with the FH;CN~ pair of complexes being more strongly 
bound. 

For both sets of isomers the complex which is hydrogen bonded through the 
N end of CN - has a much larger dipole moment. This situation occurs because, in 
all cases, the negative end of the dipole moment of the dimer is the CN - end of the 
complex and so the dipole moment of the complex is greater in magnitude when the 
negative end of the CN - moiety is farthest from the center of mass of the dimer. In 
addition, the difference between the dipole moment of FH 'NC - and FH ‘CN~ is 
0.63 D (CISD), almost exactly the CISD dipole moment of CN - (0.60 D). 

As is usual for hydrogen-bonded complexes, the dipole moment of the complex 
is greater than the vectorial sum of the two monomers. However, in this case the 
increase is much larger than normal and is partially due to the charged nature of the 
complex and the large change in the relationship between the center-of-mass and 
the center-of-electron charge which occurs upon formation of the complex. Also, 
the large polarizability of CN - probably contributes to the sizable dipole moment 
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of the complexes due to polarization of the CN electron cloud away from the HF 
or H 2 O species. 

Structures 

Table 5 lists the ab initio equilibrium structures of the FH "CN - and 
FH "NC - anionic hydrogen-bonded complexes, with those of the H 2 O 'CN and 
H 2 0 'NC - dimers presented in Table 6. See Figures 1 and 2 for the definition of 
the geometrical parameters contained in Table 6. 

For the FH"'CN - ; FH 'NC - pair of complexes the MP2 level of theory greatly 
overestimates the effects of electron correlation and CISD underestimates the impor- 
tance of electron correlation; consistent with the results obtained for the monomers. 
The CCSD method predicts bond lengths which are between the MP2 and CISD 
values, but which are much closer to the CISD values than the MP2 quantities. 
This indicates the importance of electron correlation. 

Comparing the complex tph and tcn with the monomer bond lengths we note 
that the H-F bond distance increases, as is typical upon hydrogen bond formation, 
but that the C-N bond distance decreases relative to CN - . This effect may be due 
to the loss of some of the C-N antibonding character and is supported by the earlier 
observation that the C-N distance in CN - is longer than in either HCN or HNC. 
This explanation is also consistent with the experimentally observed blue shift in 
the C-N stretch frequency of HCN "HF 54 . Since electron density is drawn away 
from the C-N bond, the C-N antibonding character is reduced leading to a shorter 
C-N distance in the complex. The C-N stretch frequency of HCN often exhibits a 
blue shift in neutral hydrogen bonded complexes 55 . 

Interestingly, although the heavy-atom distance (Rp..c or R-F.jv) I s smaller for 
FH " NC - (due to the shorter hydrogen bond), the C-N distance is more affected 
(relative to CN - ) in FH " CN - . This result tends to suggest that C contributes 
more to the C-N antibonding characteristics. Another noteworthy feature of the 
heavy-atom distances, is that for both isomers the CCSD level of theory predicts 
the largest distances while the MP2 method dramatically underestimates the heavy- 
atom lengths. 

For the H 2 0'''CN - and H 2 0"'NC - pair of complexes, CCSD geometry opti- 
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mizations were not possible. However, based upon the above comparisons between 
CCSD, CISD and MP2 for the FH;CN“ pair, it is reasonable to expect that the 
CCSD equilibrium structures will be intermediate between the CISD and MP2 opti- 
mum geometries and probably somewhat closer to the CISD structures. The various 
equilibrium structures of the H 2 0"CN- and H 2 0~NCT anionic complexes given 
in Table 6 exhibit tendencies similar to those reported above for the FH;CN~ pair 
of dimers. The O-Hj bond distance (where H a is involved in the hydrogen bond, 
see Figures 1 and 2) elongates upon complexation while the C-N linkage decreases 
relative to that in CN“. In addition, the N-Hj hydrogen bond distance is again 
shorter than the C-Hj hydrogen-bond distance; in this case by 0.138 A. The differ- 
ential heavy atom distance (i.e., R co - P -NO = 0.1402 A, CISD) is also larger than 
for the FH;CN“ pair. 

Unique to the H 2 0;CN _ complexes is the decrease in the 0-H 2 bond distance, 
the closing of angle 7 (see figures 1 and 2 for the definition of 7 ) and the non-linear 
hydrogen bond (i.e., A-H -B do not lie in a straight line). The decrease in the 0-H 2 
bond distance seems natural due to the longer 0-H a distance, though, this result 
seems to imply that the electron density of the H 2 0 monomer unit is polarized 
towards the CN - monomer unit. While this phenomenon would be expected for 
neutral hydrogen-bonded complexes, it is not necessarily expected for the case where 
one of the monomers is an anion. However, the shorter 0-H 2 distance may be 
related to the decrease in the O-H-O angle 7 . In other words, long range attractive 
forces between the H and the electron cloud around the C and/or N will result 
in a decrease in both t 0 -h 2 and 7- Such lon 8 ran S e attractions also explain the 
non-linear hydrogen bond. 

Aside from the hydrogen bond distance and the associated heavy atom dis- 
tance, the main structural difference between H 2 0"’CN - and H 2 0 "NC is the 
angle a. The smaller angle a for H 2 0 'NC - represents a larger deviation from 
linearity and is consistent with long range attractive forces between H 2 and the 
electron cloud around C (N in the case of H 2 0"CN _ ). Since the electron density 
around the C end of CN - is more diffuse, there is a stronger interaction between 
H 2 and C in H 2 0"'NC - than between H 2 and N in H 2 0"'CN . Thus, the angle a 
is smaller by about 4° for H 2 0"'NC~. 
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Diagonalization of the mass-weighted Hessian matrices explicitly demonstrate 
that each of the four complexes represents a true minimum on the PES. Therefore, 
attempts were made at the SCF level of theory to locate the transition structure 
between FH 'CN - and FH " NC - . However, due to the nature of interactions be- 
tween two closed-shell monomers the potential energy surface is very flat in this 
region in several degreees of freedom, though the total energy does rise as the CN~ 
moiety rotates. Performing the full geometry optimization is somewhat complicated 
and so the actual stationary point structure of the transition state was not pursued 
further. 


A search of the potential energy surface along the proton transfer coordinate 
was also performed in order to determine whether a second minimum (corresponding 
to F“ "HCN or F _ "HNC) exist. The search along the PES in this coordinate is 
significantly easier since all the atoms were constrained to be collinear. However, a 
second minimum (and corresponding transition state) could not be located. A brief 
discussion of the nature of the PES along the proton transfer coordinate is in order. 
Generally a transition state (TS) on a PES arises due to an avoided crossing of 
two states of the same symmetry. Thus, the SCF method is often not an adequate 
reference function for the TS. However, there are many different types of avoided 
crossings and in this particular case the SCF wavefunction should be a reasonable 
reference. This situation arises due to the fact that the orbital occupations of the 
reactants (A~"'HB) and products (AH ’B - ) are the same. What does occur as 
the proton is transferred is that two reactant molecular orbitals (a lone pair MO 
on A - and a bonding MO in HB) change character and become two product MO’s 
(a lone pair on B~ and a bonding MO in AH). However, since these MO’s are 
of the same symmetry, the transition from reactant to product MO’s is smooth 
along the proton transfer coordinate. We note that the reactant and product MO’s 
belonging to the same irreducible representation is a necessary but not sufficient 
condition for a smooth transition. Nevertheless, in this specific case the transition 
from reactant to product MO’s appears to be smooth. Therefore, the SCF function 
should represent a reasonable reference from which to evaluate dynamical electron 
correlation effects. 
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Vibrational Spectra 


The harmonic vibrational frequencies and infrared intensities for the four an- 
ionic hydrogen-bonded dimers included in this study are presented in Table 7. The 
experimental fundamentals which Ault 28 measured in matrix isolation IR studies 
are included for comparison, though, because of the rather large anharmonicities 
which the stretch modes are expected to exhibit, near quantitative accuracy with 
harmonic frequencies is not possible. The most astonishing result from Table 7 is 
the variation of the harmonic frequencies Wl FH~CN“ and Wl FH”NCT with re- 
spect to level of theory. Note that the normal mode associated with corresponds 
to the proton transfer coordinate (i.e., AH~B" - A" "HB), which for the FH;CN - 
complexes is predominately the H-F stretch. Quite clearly, an adequate treatment 
of electron correlation is extremely important in properly describing the shape of 
the potential energy surface along this coordinate. Interestingly, the large varia- 
tions in Wl (e.g., for FH" CN“, 3320 cm' 1 SCF, 2352 cm' 1 MP2 and 2844 cm' 1 
CISD) would probably not have been predicted based upon the different equilib- 
rium H-F bond distances (0.947 A SCF, 1.010 A MP2, and 0.977 A CISD), though, 
not surprisingly, there is a strong correlation between the H-F bond distance and 
the harmonic frequency. In fact, the nearly linear relationship between wjand r H F 
allows the CCSD to be estimated as ~ 2689 cm' 1 for FH CN~and ~ 2912 
cm' 1 for FH-NC". 

Given the large variation of wj with level of theory, it may seem very difficult 
to arrive at a reliable theoretical prediction for the fundamental band center vi . 
However, studies on similar systems have demonstrated that the individual har- 
monic frequency and anharmonic correction quantities converge much more slowly 
(with respect to level of theory) than does the combination, i.e., the fundamental 
band center. For example, in the study 13 of FHF by Janssen et al. the harmonic 
frequency of the antisymmetric stretch u> 3 varies from 627 cm 1 to 1538 cm 1 while 
the fundamental v$ varies only from 1427 cm -1 to 1703 cm 1 . A possible expla- 
nation for this observation may be that the A •"H" B system should be viewed as a 
particle in a one-dimensional box, where the distance R a b defines the box in which 
the proton is allowed to move. Thus, we may expect the distance Rab to converge 
more quickly (with respect to level of theory) than r ah and tbh- l n reexamining 
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the theoretical structures in Table 5 we note that once an iterative electron correla- 
tion procedure is used, then the above conditions are met (i.e., AR < Ar). In any 
case the above explanation seems feasible and will no doubt be scrutinized as more 
theoretical studies concerned with the prediction of the fundamental vibrational 
frequencies of this type of system are performed. 

The second most striking feature of the IR spectrum of the FH;CN" pair 
of complexes is the extremely large intensity exhibited by u x . Though a large 
IR intensity is expected for a mode which corresponds to proton transfer, the IR 
intensities of FH 'CN - and uq FH"NC - are even larger than the IR intensity 
reported for the analogous mode of F -- "H 2 0. However, the IR intensity reported 
for the asymmetric stretch of FHF is substantially larger than the uq FH,CN 
quantities. Interestingly, there appears to be a direct correlation between the IR 
intensity of the proton transer mode and the binding energy of the dimer. The 
appropriate IR intensity and the ab initio binding energy both decrease in the 
order FHF - > FH;CN - > F - "H 2 0 > H 2 0;CN - . This correlation suggest that 
the larger the anionic dimer binding energy then the flatter the potential energy 
surface along the proton transfer coordinate leading to a larger amplitude motion. 
Of the remaining FH;CN — normal modes, and possibly u 3 should be observable 
with u > 2 of FH-NC" also a possibility. The IR intensities of the H 2 0;CN - pair of 
complexes are more evenly distributed and, therefore, there are several vibrational 
modes which should be observable. 

For the H 2 0;CN" pair of complexes the proton transfer vibrational mode is 
u> 2 . The variation of u ? 2 with respect to level of theory is much smellier than was 
exhibited by the FH"‘CN - and FH "NC - pair, though it is still substantial. For 
example, w 2 for H 2 0"CN - is 3786 cm -1 , 3174 cm -1 , and 3497 cm -1 for the SCF, 
MP 2 and CISD levels of theory, respectively. The smaller variation of u ; 2 H 2 0;CN 
relative to uq FH;CN - was, however, to be expected due to the smaller binding 
energy of the H 2 0;CN - pair. 

Another manifestation of the smaller binding energy of the H 2 0' •CN - , 
H 2 0 “NC - pair of complexes is the lower C-N stretch harmonic frequency rela- 
tive to the FH— CN", FH"NC" pair. As noted previously, the harmonic frequency 
of CN" is less than w 2 in HCN. Thus, in an analogous manner the lower C-N stretch 
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frequency in the H 2 0;CN" pair is consistent with a smaller interaction between the 
H 2 0 and CN - monomers than exists between the HF and CN monomers. These 
observations also suggest that there is a smaller degree of charge transfer in the 
H 2 0;CN" complexes than present in the FH;CN _ pair. However, the above obser- 
vations do not indicate the relative importance of charge transfer in the bonding 
mechanism. 

Not surprisingly, the C-N stretch harmonic frequency of all four complexes 
exhibits a blue shift relative to the C-N stretch in HCN, CN" and HNC. As discussed 
earlier, neutral HCN hydrogen-bonded complexes (such as HF ’HCN) often exhibit 
a blue shift in the C-N stretch due to the loss of C-N antibonding character upon 
complexation. 

By comparing the vibrational spectra of the two FH;CN complexes or the 
H->0;CN" pair, it is evident that it would be difficult to distinguish between the 
two isomers based upon the vibrational frequencies alone. However, due to the 
differences in their structures, the best method of distinguishing the two isomers 
will be via analysis of a ro-vib rationed band. The rotational constants presented in 
Table 8 confirm this hypothesis since the differences are well within the accuracy of 
high resolution spectroscopy. 

In order to predict accurately the fundamental band centers of the vibra- 
tional modes of these complexes, a potential energy function including very high 
orders (e.g., octic terms) in some of the degrees of freedom would be required. In 
addition, a high level approach to the solution of the nuclear Schrodinger equa- 
tion, which explicitly accounts for large anharmonic couplings, would be necessary. 
This procedure would obviously be very expensive and is beyond the scope of the 
present study. However, the vibrational analysis that we have performed has lead 
to further insight concerning the proton transfer vibrational mode, the most likely 
fundamental of the FH;CN~ and H 2 0;CN* anionic dimers to be experimentally 
observed. Furthermore, the similarity of the vibrational spectra of the pairs of iso- 
mers has been explicitly demonstrated and a method by which the isomers may be 
spectroscopically distinguished has been noted. 


Bonding 
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The much larger binding energies found in anionic hydrogen-bonded com- 
plexes (relative to neutral hydrogen-bonded complexes) lead to questions concerning 
the nature of this interaction. For example, if Morokuma and coworkers’ hydrogen 
bond energy decomposition scheme is applied, which components exhibit signifi- 
cantly different characteristics for anionic complexes? As discussed earlier, such an 
analysis has been performed 17 on FHF - , however, it seems likely that an asym- 
metric anionic hydrogen-bonded complex will possess quite different characteristics 
than FHF - where charge transfer is clearly very important. Moreover, the binding 
energy of FHF - (~ 39 kcal/mole) is significantly larger than that for the complexes 
included in this study. 

The three hydrogen bond components which one might intuitively expect to 
yield large attractive energies are the electrostatic, polarization and charge trans- 
fer interactions, i We will not discuss the electrostatic interactions here except 
to note that the detailed structure of this interaction must be very different for 
FH'"CN~ (H 2 0 ’CN - ) and FH’NC - (H 2 0"NC - ) because of the reversal of the 
dipole moment of CN - . The total binding energies are very similar, however. Even 
though the dipole of CN - has been reversed, this does not mean that the toted 
electrostatic energies of the two isomers are different, though it does seem probable 
that there will be a detectable difference. In the latter case, some other hydrogen 
bond energy component must compensate. 

The polarization interaction for the anionic hydrogen-bonded complexes in- 
cluded in this study must be significantly larger than exists in a similar neutral 
hydrogen-bonded complex. This conclusion is based upon the much larger polar- 
izability which anions possess (e.g., at the SCF TZ2P+diffuse level of theory the 
mean polarizability for HF, H 2 0 and CN - is 0.66, 1.14 and 3.48 A 3 , respectively). 
Furthermore, the decomposition analyses which have been performed on neutral 4 
and anionic 17 (FHF - ) hydrogen-bonded complexes provide additional support for 
this inference. 

It is difficult to assess the degree of charge transfer. One method would be 
to perform a Mulliken population analysis on the complex, and from these data 
determine the number of electrons associated with each monomer. Performing such 
an analysis on the FH 'CN - and FH 'NC - anionic complexes and comparing to 
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a similar analysis on the HCN ' HF hydrogen bonded complex shows that indeed 
there is more charge transfer in the anionic species. However, as is well known, a 
Mulliken population analysis associates electrons to a given nucleus in an ambiguous 
manner. Therefore, given the rather small differences between the neutral and 
anionic complexes the validity of the results would seem to be in question. An 
alternative method would be to perform electron density difference plots between 
the complexes and their respective monomers and compare these for the anionic 

an d neutral hydrogen-bonded complexes. 

Valence electron density difference plots from CISD natural orbitals have been 
performed for the FH'CN - , FH~NC“ and HCN’HF hydrogen-bonded complexes 
and are presented in figures 3 - 5, respectively. The contour interval for all three 
plots is the same. The HCN ' HF equilibrium geometry was taken from reference 
9, but the TZ2P + diffuse basis set of the current study was used. Noting that 
short dashed lines indicate electron depletion and solid lines indicate an increase 
in electron density, it is clear that the anionic complexes exhibit a larger charge 
transfer from the CN~ species to the HF monomer than occurs in the neutral 
complex. Moreover, this conclusion is enforced by the large buildup of electron 
density behind the F atom in the anionic complexes. Interestingly, the plots also 
show a depletion of electron density in the C-N bonding region upon complexation. 
This observation is entirely consistent with earlier statements concerning the C-N 
equilibrium bond length and harmonic frequency in HCN, HNC and CN . 

Considering the above discussion and previous results 2 4,17 , a possible scene- 
rio may be suggested. It is likely that the electrostatic, polarization and charge 
transfer energy components of an anionic complex are all larger than those for a 
similar neutral complex. Moreover, as one progesses from an asymmetric complex to 
a symmetric species (i.e., the proton half way between the heavy atoms), the charge 
transfer component will become much larger. This model, then, also explains the 
large difference between the binding energies of FHF" and FH;CN _ . 

Concluding Remarks 

The FH "CN _ and FH 'NC - pair of anionic hydrogen-bonded complexes have 
been shown to be nearly isoenergetic and the theoretical binding energy is in good 
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agreement with experiment. The H 2 O 'CN“ and H20"'NC” pair of complexes 
are also very close energetically with the best ab initio binding energy again in 
good agreement with the experimental value. The equilibrium structures of the 
isomers, however, do exhibit small differences (e.g., the N "H hydrogen bond is 
shorter than the C"H hydrogen bond) which lead to slightly different rotational 
constants. Thus, because the harmonic IR spectra of the two pairs of isomers are 
so similar, the different rotational constants provide a means by which the isomers 
may be experimentally distinguished. It is concluded, however, that an accurate 
theoretical determination of the fundamental frequencies will require a large portion 
of the potential energy surface to be investigated using a high level of electronic 
structure theory, such as CCSD coupled with a large one-particle basis set. In 
addition, a sophisticated solution of the nuclear motion problem capable of treating 
large anharmoniGities will be necessary. 

Another significant outcome of this study involves the CCSD investigations of 
the monomers. This is the first study which has fully optimized molecular structures 
and evaluated several equilibrium molecular properties at the CCSD level of theory 
with a large one-particle basis set (i.e., larger than double zeta plus polarization) 
for chemical systems exhibiting a range of bonding characteristics. The CCSD 
equilibrium structures, harmonic frequencies, dipole moments and IR intensities for 
HF and H 2 0 clearly demonstrate that near quantitative results may be obtained 
for systems which are well described by a single determinant reference function. 
Although the CCSD results for HCN, HNC, and OH" have slightly larger errors, 
they are still very good and are superior to the analogous MP2 and CISD quantities. 
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Figure Captions 

Figure 1. Definition of the internal coordinates for the H20 "CN - anionic hydrogen- 
bonded complex. 

Figure 2. Definition of the internal coordinates for the H20‘"NC~ anionic hydrogen- 
bonded complex. 

Figure 3. Valence electron density difference plot for the FH'"CN“ anionic complex. 
Short dashed lines indicate a depletion of electron density while solid lines indicate 
an increase of electron density. 

Figure 4. Valence electron density difference plot for the FH‘”NC - anionic complex. 
Short dashed lines indicate a depletion of electron density while solid lines indicate 
an increase of electron density. 

Figure 3. Valence electron density difference plot for the HCN'HF hydrogen- 
bonded complex. Short dashed lines indicate a depletion of electron density while 
solid lines indicate an increase of electron density. 
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Table 1 

Theoretical predictions of the total energy, optimum structure 
and dipole moment of the possible fragmentation monomers. 
Energies, bond lengths, angles and dipole moments are given in 
Hartrees, A, degrees and Debyes respectively. 


Monomer 

Structure 

Method 

Energy 

Dipole Moment 

F“ 



SCF 

-99.455226 

- 




MP2 

-99.724127 

- 




CISD 

-99.691431 

- 




CCSD 

-99.718898 

- 

HF 

*HF 

0.8985 

SCF 

-100.064599 

1.89 



0.9208 

MP2 

-100.317487 

- 



0.9149 

CISD 

-100.293761 

1.82 



0.9176 

CCSD 

-100.319608 

1.81 



0.9168 

Expt° 

- 

1.80 

OH" 

*OH 

0.9426 

SCF 

-75.413857 

1.35 



0.9655 

MP2 

-75.688586 

- 



0.9580 

CISD 

-75.654436 

1.25 



0.9624 

CCSD 

-75.685247 

1.24 



0.9643 

Expt 6 

- 

" 

CN- 

*cn 

1.1513 

SCF 

-92.342662 

0.46 



1.1870 

MP2 

-92.695252 

- 



1.1688 

CISD 

-92.637209 

0.60 



1.1744 

CCSD 

-92.698406 

0.64 



- 

Expt 

- 

- 

h 2 o 

ro/f 

0.9404 

SCF 

-76.062199 

1.95 


ZHOH 

106.4° 





*OH 

0.9593 

MP2 

-76.315821 

- 


Z.HOH 

104.5° 





toh 

0.9540 

CISD 

-76.293736 

1.89 


ZHOH 

105.0° 





row 

0.9571 

CCSD 

-76.321837 

1.88 


ZHOH 

104.8° 





toh 

0.9578 

Expt c 

- 

1.85 


ZHOH 

104.5° 

- 




Table 1 continued 


Monomer 

Structure 

Method 

Energy 

Dipole Moment 

HCN 


1.0572 

SCF 

-92.909715 

3.27 


tatc 

1.1235 





THC 

1.0641 

MP2 

-93.262244 

- 


r nc 

1.1637 





*HC 

1.0618 

CISD 

-93.207319 

3.08 


tNC 

1.1438 






1.0653 

CCSD 

-93.268033 

3.03 


rjvc 

1.1502 





*HC 

1.065 

Expt d 

- 

2.99 


*NC 

1.153 




HNC 

r hn 

0.9819 

SCF 

-92.892533 

2.97 


r cn 

1.1440 





r hn 

0.9959 

MP2 

-93.233324 

- 


*CN 

1.1733 





r hn 

0.9903 

CISD 

-93.185062 

3.09 


r ex 

1.1609 





IHN 

0.9939 

CCSD 

-93.244549 

3.11 


*CN 

1.1658 






0.9940 

Expt e 

- 

3.05 


*CN 

1.1689 



- 


a. All experimental structures refer to derived equilibrium structures. HF bond length from Ref [56] and 
dipole moment from Ref [57]. 

b. Ref [58]. 

c. Structure from Ref [59] and dipole moment from Ref [60]. 

d. Structure from Ref [61] and dipole moment from Itef [62]. 

e. Structure from Ref [63] and dipole moment from Ref [64]. 


Table 2 

Harmonic vibrational frequencies and infrared intensities 
for the monomers. Frequencies are given in cm" 1 and IR intensities 
(in parentheses) in km/moL 


Monomer 

Normal 

mode 

SCF 

MP2 

CISD 

CCSD 

Expt 

HF 

Wl(ff) 

4469 (164) 

4126 

4212 (114) 

4165 (106) 

4139 (96)° 

CN~ 

Wi(o-) 

2317 (45) 

1982 

2167 (21) 

2112 (16) 

- 

OH" 

Wi(<t) 

4073 (62) 

3805 

3855 (77) 

3782 (85) 

3738* 

H 2 0 

wi(ai) 

4130 (15) 

3841 

3919 (6) 

3865 (4) 

3832 (2) e 


w 2 (ai) 

1757 (96) 

1657 

1694 (74) 

1684 (71) 

1649 (54) 


w 3 (b 2 ) 

4233 (92) 

3967 

4021 (62) 

3972 (57) 

3942 (45) 

HCN 

w lW 

3608 (72) 

3451 

3497 (68) 

3438 (64) 

3442 (59) d 


w 2 (<r) 

2407 (11) 

2027 

2236 (2) 

2171 (0.4) 

2129 (0.2) 


w 3 (tt) 

855 (70) 

686 

734 (72) 

706 (72) 

727 (50) 

HNC 

Wl((T) 

4046 (379) 

3818 

3899 (286) 

3839 (260) 

3842* 


w 2 (<t) 

2282 (103) 

2017 

2145 (70) 

2098 (62) 

2067 


w 3 (tt) 

472 (313) 

459 

425 (277) 

442 (269) 

490 


a. All experimental frequencies are derived harmonic frequencies. The harmonic frequency is taken from ref 
[56] while the IR intensity is taken from ref [65]. 

b. The harmonic frequency is taken from ref [58]. 

c. The harmonic frequencies are taken from ref [66] while the IR intensities are taken from ref [67]. 

d. The harmonic frequencies are taken from ref [68] while the IR intensities are taken from ref [69]. 

e. The harmonic frequencies are taken from ref [63]. 


Table 3 

Thermochemical data for possible fragmentation products of the 
titled anionic hydrogen bonded complexes in kcal/mole. 


Method 

AH 298 “- c 

AH 298 ft 

SCF 

27 

50 

MP2 

17 

36 

CISD 

21 

42 

CCSD d 

20 

41 

Expt e 

26 

45 


a) HF + CN~ + Aif 2 98 -*■ F ~ + HCN 

b) H 2 0 + CN~ + AH 2 9» -> 0H~ + HCN 

c) Difference in total electronic energies, zero point vibrational energies and rotational and 
translational contributions at 298K. 

d) The CCSD values use the CISD zero point vibrational energies. 

e) Ref [70]. 


Table 4 

Predicted binding energies (kcal/mol) and dipole moments (Debyes) 
for the anionic hydrogen bonded complexes. The binding energies 
were computed with respect to the most stable dissociation products 
as indicated in Table 3. 


Anionic 

complex 

Method 

Energy 

AE a 

AE h 

Dipole 

Moment 

FH CN" 

SCF 

-192.442703 

22.1 

22.0 

2.92 


MP2 

-193.056400 

27.8 

26.8 



CISD 

-192.940474 

25.3 

24.3 

2.38 


CCSD 4 

-193.059361 

25.5 

24.5 

- 

FH-NC" 

SCF 

-192.444070 

23.0 

22.9 

3.19 


MP2 

-19.3.054960 

26.1 

25.0 

- 


CISD 

-192.940721 

25.4 

24.3 

3.01 


CCSD 4 

-193.059157 

25.4 

24.3 

- 

- 

Expt d 

- 

- 

21.1 

“ 

H 2 0-CN- 

SCF 

-168.424382 

12.4 

12.3 

4.29 


MP2 

-169.036450 

16.3 

15.6 

- 


CISD 

-168.920387 

14.5 

13.8 

3.84 


CCSD* 

-169.043811 

15.2 

14.5 


H 2 0 NC- 

SCF 

-168.425723 

13.4 

13.4 

4.48 


MP2 

-169.036487 

16.2 

15.4 

- 


CISD 

-168.921256 

15.1 

14.3 

4.32 


CCSD 4 

-169.044430 

15.5 

14.7 

- 


Expt-^ 



12.7±0.8 



a. Includes zero point energy and translational, rotational correction for 298K, see ref 71 for method. 

b. Includes zero point energy and translational, rotational correction and basis set superposition error 

determined by the counterpoise method, ref 53. .. . 

c. CCSD energy performed at CCSD equilibrium geometry. The single point energy allowed all or 1 a 
to be active, whereas in the geometry optimization the core and corresponding virtual orbitals were frozen. 
The optimum CCSD energies are -193.004605 and -193.004374 for FH CN* and FH NC~, respectively. 
CISD zero point energies were used. 

d. Ref 24. . , ^cn 

e. CCSD energy at the CISD equilibrium geometry. All orbitals active m the CCSD procedure. CISD zero 

point energies were used. 

f. Ref 29. 


Table 5 

Goomotric.nl structures for the anionic hydrogen bonded complexes 
FH "CN - and FH'NC - . Bond lengths are given in A. 


Anionic complex 

Method 

r fii 

TCN 


R-F...C 

FH-CN- 

SCF 

0.9472 

1.1457 

1.7790 

2.7262 


MP2 

1.0103 

1.1803 

1.6039 

2.6142 


CISD 

0.9769 

1.1585 

1.6768 

2.6537 


CCSD 

0.9879 

1.1696 

1.6666 

2.6546 

- 


rpu 

*CN 

rir...N 

R -F...N 

FH-NC- 

SCF 

0.9437 

1.1483 

1.6481 

2.5918 


MP2 

0.9898 

1.1820 

1.5460 

2.5358 


CISD 

0.9678 

1.1611 

1.5798 

2.5476 


CCSD 0.9768 1.1719 1.5785 2.5553 



Table 6 

Geometrical structures for the anionic hydrogen bonded complexes 
H 2 0 CN” and H 2 0 NC~. See Figure 2 for definitions of 
the molecular bond angles. Bond lengths are given in A 
and bond angles in degrees. 


Anionic 

complex 

Method 

*OH x 

r oh-2 

*CN 

rc.^i 

R CO 

Ot 

0 

7 

H 2 0 CN" 

SCF 

0.9599 

0.9392 

1.1484 

2.1246 

3.0706 

176.3° 

168.2° 

103.6° 


MP2 

0.9966 

0.9586 

1.1836 

1.9055 

2.8984 

t— ‘ 
-4 

bo 

o 

173.9° 

102.0° 


CISD 

0.9769 

0.9486 

1.1614 

1.9903 

2.9601 

176.6° 

171.5° 

102.7° 



* oh x 

r oh 7 

r cn 


R/vo 

a 

0 

7 

HjO'NC- 

' SCF 

0.9596 

0.9390 

1.1500 

1.9501 

2.8963 

173.0° 

168.3° 

103.6° 


MP2 

0.9914 

0.9581 

1.1847 

1.8012 

2.7885 

172.8° 

173.6° 

102.2° 


CISD 0.9749 0.9484 1.1629 1.8523 2.8199 172.3° 171.4° 102.8 ( 


Table 7 

Harmonic vibrational frequencies and infrared intensities 
for the anionic hydrogen bonded complexes. 

Frequencies arc given in cm -1 and IR intensities (in parentheses) in km/mole . 


Anionic 

Complex 

Normal 

mode 

SCF 

MP2 

CISD 

Expl,.* 

FH- 

•CN- 


3320 (2203) 

2352 

2844 ( 2679) 

1800 



uj 2 (( 7 ) 

2373 

(21) 

2044 

2261 

(7) 

2500 



t) 

246 

(46) 

298 

278 

(61) 

- 



W 4 (7t) 

1043 

(190) 

1142 

1107 

(143) 

1100 



w s (tt) 

162 

(3) 

159 

163 

(3) 

- 

FH- 

■NC" 


3406 ( 

(2271) 

2727 

3034 (2618) 

1800 



u 2 (<r) 

2346 

(101) 

2030 

2235 

(74) 

2500 




276 

(58) 

309 

300 

(70) 

- 




1034 

(254) 

1104 

1086 

(211) 

1100 



w s (?r) 

123 

(5) 

127 

126 

(8) 

- 

h 2 o- 

•CN- 

wi(a/) 

4195 

(43) 

3911 

4050 

(28) 





3786 

(801) 

3174 

3497 

(1106) 




a>s(a/) 

2345 

(30) 

2012 

2232 

(12) 




o.> 4 (a/) 

1822 

(114) 

1722 

1777 

(88) 




ws(a/) 

406 

(73) 

472 

446 

(66) 




Wfi(a/) 

175 

(26) 

217 

202 

(32) 




w 7 (a/) 

96 

(^) 

93 

98 

(3) 




w«(a//) 

798 

(116) 

902 

860 

(84) 




L>g (n//) 

115 

(13) 

106 

114 

(10) 


h 2 o 

-NC- 

wi(a/) 

4198 

(43) 

3919 

4054 

(30) 




<x> 2 ( a/) 

3800 

(864) 

3289 

3548 i 

(1116) 




w,3(a/) 

2329 

(74) 

2004 

2216 

(50) 




w 4 (a/) 

1827 

(127) 

1730 

1782 

(100) 




cj 5 (n/) 

398 

(73) 

463 

439 

(66) 




w fi (a/) 

197 

(34) 

237 

225 

(41) 




w 7 (a/) 

71 

(4) 

75 

72 

(5) 




Vfi(a//) 

802 

(133) 

890 

856 

(101) 




UJg(a//) 

90 

(30) 

96 

93 

(33) 



a. Experimental fundamental frequencies are taken from ref [28]. 


Table 8 

Rotational Constants (MHz) for the Equilibrium structures of the 

Anionic Complexes. 


Anionic SCF MP2 CISD CCSD 

Complex 


FH-CN" 

A 

3825 

4021 

3966 

3946 

FH-NC" 

A 

H 2 0-CN- 

4348 

4438 

4443 

4404 

' A 

513050 

536470 

530720 


B 

3370 

3631 

3548 


C 

H 2 0-NC- 

3348 

3607 

3524 



A 

B 

C 


455610 

3887 

3854 


469440 

4065 

4031 


455490 

4033 

3998 


Figure Captions 


Figure 1. Definition of the internal coordinates for the H20"'CN anionic hydrogen- 
bonded complex. 

Figure 2. Definition of the internal coordinates for the H 2 O "NC- anionic hydrogen- 
bonded complex. 

Figure 3. Valence electron density difference plot for the FH'"CN~ anionic complex. 
Short dashed lines indicate a depletion of electron density while solid fines indicate 
an increase of electron density. 

Figure 4. Valence electron density difference plot for the FH NC - anionic complex. 
Short dashed fines indicate a depletion of electron density while solid fines indicate 
an increase of electron density. 

Figure 5. Valence electron density difference plot for the HCN 'HF hydrogen- 
bonded complex. Short dashed fines indicate a depletion of electron density while 
solid fines indicate an increase of electron density. 
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